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Mediante simulaciones de Monte Carlo, utilizando escalamiento de tamano finito y transforma- 
ciones conforme se reportan evidencias de una transicion de fase topologica en cristales liquidos en 
dos dimensiones. A altas temperaturas se presenta una fase paramagnetica mientras que a bajas 
temperaturas se encuentra un fase de cuasi-largo-orden (QLRO). Se encuentra que a muy bajas tem- 
peraturas el exponente de la funcion de correlacion de la fase QLRO es lineal con la temperatura 
teniendose de esta manera un comportamiento tipico de ondas de espin. Esto ultimo contradice 
predicciones de que para sistemas con grupo de simetria global no abeliano las ondas de espin no 
son relevantes. Adicionalmente se discute que implicaciones tendria la presencia de impurezas en 
este tipo de transiciones. 
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INTRODUCCION 

Mermin y WagnerQ establecieron que para sistemas con grupo de simetria continuo no existe fase ferromagnetica 
o de orden a largo alcance (LRO), a temperatura diferente de cero, para dimensiones menores o iguales que dos. 
Sin embargo, sistemas de este tipo pueden presentar un tipo de transicion que esta determinado por la aparicion 
de defectos topologicos en pares a bajas temperaturas que justo en la transicion, Tkt, se desliganj^, HI Q ■ A este 
fenomeno se le llama transicion de fase topologica o transicion de Berezinskii, Kosterlitz y Thouless (BKT). 

El modelo XY, con grupo de simetria 0(2), presenta este tipo transicion en dos dimensiones 4]. Entrc sus carac- 
teristicas resaltantes tenemos que la fase LRO es sustituida, a bajas temperaturas, por un orden a cuasi-largo-alcance 
(QLRO). A muy bajas temperaturas, las correlaciones son dominadas por las ondas de espin obteniendose una de- 
pendencia del exponente de la funcion de correlacion, 77, con la temperatura, T, de la forma: r/ = T/2n, dondc 
la constante de Boltzmann ks — 1 y el factor de acoplamiento J = 1. Debido a este comportamiento QLRO, a 
temperaturas menores que Tkt, l a susceptibilidad, que mide las fluctuaciones de la magnetization, diverge. Otra 
caracteristica de este tipo de transicion es que a temperaturas justo por encima de la temperatura de transition 
t = (T — Tkt)/Tkt ft 0, la longitud de correlacion, £, diverge de una manera mucho mas fuerte que la tipica ley de 
potencias, £ ~ t~ u , encontrada en las transiciones de segundo orden. Dicha divergencia es del tipo de una singularidad 
escencial: £ ~ exp(&t -1 / 2 ). 

Por otro lado, otro modelo de simetria continua como el de Heisenberg ferromagnetico, con grupo de simetria 0(3), 
no presenta ningun tipo de transicion de fase en d = 2 mientras que si presenta la tipica transicion para-ferromagnetica 
en d = 30,0]. Por esta razon surgio la pregunta de que si sistemas con grupo de simetria no abelianos podian presentar 
transiciones del tipo BKT. Se ha reportado que el modelo de Heisenberg antiferromagnetico completamente frustrado 
presenta una transicion del tipo BKT pero sin correlaciones del tipo de ondas de espin a bajas temperaturas 0,0] • 

Kunz y ZumbachJ^ realizaron un estudio intensivo en d = 2, mediante simulaciones, del modelo RP 2 el cual consiste 
en un sistema con grupo de simetria global 0(3) pero con grupo de simetria local Z-i- El modelo RP 2 describe la 
transicion nematica isotropica de cistales liquidos en d — 3. Ellos determinaron la longitud de correlacion para t > y 
encontraron un buen ajuste para una singularidad escencial. Por otro lado, mediante escalamiento de tamaho finito, 
cstimaron que tambien era valido un ajuste del tipo ley de potencia. Sin embargo, en base a calculos de energia y 
calor especifico asi como de ciertas cantidades que estiman el mimero de defectos topologicos, argumentaron que la 
transicion de fase deberia ser del tipo BKT. 

En el 2003 se retomo el estudio pero ahora utilizando el modelo de Lebwohl-Lasher(LL) [Toj para describir a los 
cristales liquidos en d = 2. Este modelo ha sido muy exitoso para detectar la transicion de fase discontinua debil que 
se observa en los experimentos de cristales liquidos en d — 3. En este modelo se representan las moleculas mediante 
vectores unitarios a w colocados en los sitios de una red hipercubica A de longitud L. Para este sistema el Hamiltoniano 
viene dado por: 

H J 

-]^T = 1^tY,Y, P2 (^-^)> (!) 
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donde P2 es el segundo polinomio de Legendre y la interaction es de primeros vecinos. En d — 2, se determinojllj la 
presencia de un orden del tipo QLRO a bajas temperaturas con soluciones del tipo de ondas de espi'n cuando T tiende 
a cero. Esto se hizo mediante el novedoso metodo de las transformaciones conformes. De esta manera se concluyo, 
que sistemas con grupo de simetria no abeliano tenian transition BKT con correlaciones del tipo de ondas de espi'n a 
bajas temperaturas. 

En el presente articulo revisara el metodo de transformaciones conformes y su utilidad para determinar exponcntcs 
de la funcion correlation para sistemas con invarianza de escalas. Luego estudiando el parametro de orden nematico 
bajo Tkt-, mediante escalamiento de tamano finito, se obtendra de nuevo el exponente de la correlation y se mostrara 
la excelente concordancia entre los dos procedimientos. Finalmente se discutira el efecto de la presencia de impurezas 
sobrc cstc tipo do transition. 

METODO DE LAS TRANSFORMACIONES CONFORMES 

Uno de los grandes problemas que existe al simular un sistema fisico es el de tener evidencia de lo que ocurre 
en el lfmite termodinamico. Por lo general, en los estudios de fenomenos criticos, se simulan sistemas a diferentes 
tamanos y se estudian como escalan las cantidades termodinamicas en funcion de L. El exponente de escalamiento 
de dichas cantidades tiene relation con los exponentes criticos asociados a dichas cantidades. El costo computational 
es demasiado alto. Recientemente, para sistemas en d = 2 se ha comenzado a utilizar las transformaciones conformes 
(TC). Esta tecnica consiste en realizar simulaciones en sistemas flnitos y conectar los resultados con los de un sistema 
infinito via una TC. Esto lo podemos realizar en sistemas que presenten invarianza de escalas y para d = 2. 

Sistemas con invarianza de escala cumplen con la hipotesis de homogeneidad para la funcion de correlation de dos 
puntos para cualquier densidad 0, tal como el parametro de orden, la energia, etc., (<p(br2)<p(bri)) = {4>(r2)4>(fi)) , 
donde b es un factor de escala y rj es la dimension anomala o exponente de la correlation. De la misma manera, si existe 
invarianza de escala, podemos relacionar mediante la hipotesis de homogeneidad la funcion de correlation entre sis- 
temas que se conectan bajo una transformation conforme. Supongamos que las simulaciones se realizan en un sistema 
cuadrado (w — u + iv) de tamaho L x L (-L/2 < u < L/2, < v < L) y mediante una transformation de Schwarz- 
Christoffel lo mapeamos en el piano semi-infinito (z = x + iy, < y < 00). Para esta transformation el cambio de es- 
cala es local y las funciones de correlation se relacionan mediante ((j>(wi)<j>(w2)) = \w' ' (zi)\~ x ' t '\w l '{z2)\~ x ' 1 ' {<j)(zi)4>(z2)) , 
donde x^ = \r\. 

Si las simulaciones se realizan en una red cuadrada con condiciones de borde fijas se obtiene que el perfil de densidad 
(funcion de correlation de un punto) se comporta como una ley de potencia de la forma: (0(w)) sq ~ k(w)~ 1) / 2 con 
k(w) — $s[z(w)](\l — z(w) 2 \\l — k 2 z{w) 2 \)~ 1 ^ 2 , donde z(w) = sn(2Kw/L), siendo K la integral eliptica completa de 
primera clase, sn el seno eliptico de Jacobi y k una constante. Esta expresion se obtiene ya que el perfil para una red 
semi-infinita con condiciones de borde fijas se conoce exactamente ~ t/ - ^^ [T3| . 

Gracias a la invarianza de escalas presente en las fases del tipo qlr6, la metodologia de las TC se aplico con 
exito en la determination de los exponentes de la funcion de correlation del modelo XY en d = 2 a 
Mas recientemente, utilizando las TC en el modelo LL para los cristales liquidos en d = 2, se concluyo la existencia 
de un orden QLRO en este sistema L Lj . Para este problema la funcion de correlation de un punto utilizada es: 
m(w) = {P2{&w ■ hg\f w ^)) S q donde h,Qu w ) mdica que todas las moleculas en el borde dk(w) tienen la orientation fija 
h. En la Fig. 1(a) se muestran los perfiles del parametro de orden como funcion de k. Se nota, como a temperaturas 
menores que Tkt, se presenta un comportamiento del tipo ley de potencias mientras que muy por encima de esta 
temperatura dicho comportamiento es abandonado indicando que el sistema se encuentra en una region donde no 
existe invarianza de escalas. La determination del cambio de comportamiento se puede hacer de manera cuantitativa 
calculando el X 2 E3- En la Fig. 1(b) se grafican los exponentes de la correlation obtenidos a partir de los ajustes de 
la Fig. 1(a) como funcion de T. A muy bajas temperaturas se observa el comportamiento lineal tipico de las ondas 
de espi'n. 

ESCALAMIENTO DE TAMANO FINITO 

En el presente trabajo se realizaron simulaciones del modelo LL en dos dimensiones con condiciones de borde 
periodicas. Las simulaciones se hicieron para temperaturas comprendidas entre T — 0.4 y T = 0.58. Se simularon 
sistemas con longitudes comprendidas entre L = 16 y L = 512. Se empleo el algoritmo de Wolff0. Se tomaron del 
orden de 10 5 pasos de equilibration y 10 6 pasos de promediacion. Los tiempos de relajacion para todos los tamanos 
y longitudes no sobrepasaron los 200 pasos de Monte Carlo. En particular, se reportan estimados del exponente de la 
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Figure 1: a) Parametro de orden m(w) local en funcion de k(w) para diferentes temperaturas. Se nota un comportamiento 
del tipo ley de potencias para valores de temperatura menores que el valor Tkt = 0.513 obtenido por Kunz y Zumbach|i|- b) 
Exponente de la funcion de correlation como funcion de la temperatura utilizando los ajustes de la Fig. 1(b) junto con los 
obtenidos para L — 48 y L — 200. Los triangulos representan los estimados hechos a partir de escalamiento de tamano finite 
A temperaturas cercanas a cero se obseva el comportamiento r\ cx T. 



funcion de correlation del parametro de orden utilizando la tecnica de escalamiento de tamano finito para el parametro 
de orden. 

Un buen parametro de orden para describir la transition de fase nematica-isotropica, en d = 3, para el modclo LL 
es: 



M = L- 2 r^P 2 (a u 



■h))=L- 2 



^P 2 (cos0 w ) 




(2) 



donde h es un vector unitario que indica la direction preferential y se le denomina director y 9 W es el angulo entre a w 
y el director. M tiende a 1 a muy bajas temperaturas y a a altas temperaturas. En d = 2, M deberia ser para 
todo T ya que no puede haber rompimiento de la simetrfa continuaQ. 

En la Fig 2. (a) se tiene una grafica del parametro como funcion de la temperatura. Se observa un valor finito del 
parametro de orden a temperaturas bajas. Sin embargo, la tendencia es que disminuya al aumentar L. En la Fig. 2(b) 
se muestra el comportamiento del parametro de orden con la longitud del sistema. Observese que para temperaturas 
bajas el comportamiento es del tipo ley de potencia. Esto indica la presencia de invarianza de escalas para T < Tkt 
o lo que es lo mismo un orden del tipo QLRO. Sobre Tkt ia caida es exponential, comportamiento caracteristico de 
una fase paramagnetica. El exponente de la ley de potencia obtenido a bajas temperauras es el correspondiente a la 
funcion de correlation de un punto, es decir, ^rj. Colocando los exponentes provenientes de los ajustes realizados con 
los datos de la Fig. 2(b) en la Fig. 1(b) tenemos que la correspondencia entre las dos metodologfas reportadas es 
excelente. 



DISCUSION 



En el presente trabajo se reportan evidencias de que a bajas temperaturas existe un orden del tipo QLRO. La 
manera mas eficiente de realizar estos estimados es mediante el uso de la tecnica de las transformaciones conformes 
ya que solo se necesita hacer simulaciones en un solo tamano de red, inclusive para los primeros resultados basto el 
uso de una red pequena [L = 48). El costo computacional del metodo de escalamiento de tamano finito es mucho 
mayor al tener que simular para muchas longitudes de red y tamanos mayores. Realizando estudios del escalamiento 
de la susceptibilidad se tiene un porcedimiento alternativo para el calculo del exponente rj al igual de una forma muy 
precisa para determinar que en Tkt la singularidad es escencial0|. Por lo tanto, defmitivamente se concluye que, 
para los cristales liquidos en d — 2, se tiene un sistema con grupo de simetrfa global no-abeliano que presenta una 
transition BKT con orden QLRO a bajas temperaturas. 
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Figure 2: a) Parametro de orden en funcion de la temperatura para diferentes tamano L del sistema. Se observa como al 
incrementar L el valor de M para cada T disminuye. El valor de Tkt mostrado es el estimado por Kunz y Zumbach b) 
Para cada T el parametro de orden en funcion de L. Comportamiento del tipo ley de potencia para T < Tkt- Por encima de 
la temperatura crftica la tendencia es la de un comportamiento exponencial. 



En la literatura se ha reporta do q ue para sistemas con transicon de tipo BKT, la introduction de desorden en los 
enlaces es totalmente irrelevante 17]. Esto se debe a que debido a la singularidad escencial para t > el exponente 
del calor espetifico a — > —oo (2 — dv = a) y debido al criterio de Harris .1§| los exponentes cn'ticos del sistema puro no 
deben'an cambiar. Sin embargo, debido a cambios que aparecen en la coordination el valor de Tkt deberfa depender 
de la intensidad, c, del desorden. Si se realizaran simulaciones del LL con enlaces aleatorios obtendrfamos que el 
exponente de la correlation seria identico para todos los valores Tkt{c). 

Donde resulta muy interesante discutir el problema del desorden seria en d = 3 con campo magnetico aleatorio. 
Segun el criterio de Imry y Ma |Hj | la fase nematica deberfa desaparecer y podria ser sustituida por un orden del tipo 
QLROlIllHl. 
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